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ABSTRACT 

V440 Per is a Population I Cepheid with the period of 7.57 day and low amplitude, 
almost sinusoidal light and radial velocity curves. With no reliable data on the 1st 
harmonic, its pulsation mode identification remained controversial. We obtained a 
radial velocity curve of V440 Per with our new high precision and high throughput 
Poznah Spectroscopic Telescope. Our data reach the accuracy of 130 m/s per individual 
measurement and yield a secure detection of the 1st harmonic with the amplitude of 
A2 = 140 ± 15 m/s. The velocity Fourier phase 02i of V440 Per is inconsistent at 
the 7.25(7 level with those of the fundamental mode Cep heids, implying that the star 
must be an overtone Cepheid, as originally proposed bv lKienzle et al.l (|l999f ). Thus, 
V440 Per becomes the longest period Cepheid with the securely established overtone 
pulsations. We show, that the convective nonlinear pulsation hydrocode can reproduce 
the Fourier parameters of V440 Per very well. Requirement to match the observed 
properties of V440 Per constrains free parameters of the dynamical convection model 
used in the pulsation calculations, in particular the radiative losses parameter. 

Key virords: stars: variables: Cepheids - stars: oscillations - stars: individual: 
V440 Per - techniques: spectroscopic - methods: data analysis - hydrodynamics 



1 INTRODUCTION 

. The Population I sinusoidal or s-Cepheids are a small group 
' of the Cepheids pulsating in the first radial overtone. In the 
Galaxy, where individual Cepheid distances are usually not 
accurately known, the s-Cepheids are discriminated from 
the fundamental mode pulsat ors with the Fourier d ecom- 
position of their light curves l|Antonello et al.lll990l ). The 
method works well only for the variables with periods be- 
low 5 days. Fortunately, it can be extended to longer pe- 
riods with the help of t he ra dial velocity curves, as was 
shown by iKienzle et al.l l|l999l ). Studies of the longer pe- 
riod overtone Cepheids are of the great interest. The ve- 
locity Fourier parameters of the s-Cepheids display a very 
characteristic progression with period, which is attributed 
to the 2:1 resonance at P = 4.2 — 4.6 day between the first 
and f ourth overtones (|Kienzle et al.|[l999l : iFeuchtinger et al] 
l200d ). Unfortunately, any detailed modeling of this progres- 
sion and pinpointing of the resonance position is hampered 
by a scarcity of the s-Cepheids with P > 5.5 day, with 
MY Pup being the only secure identification. One more can- 

* Based on observations from Poznari Spectroscopic Telescope 
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dida te, V440 Per (P = 7.57 day), has been identified by 
Kie nzle et al. 

I ([1999). They noted, that the velocity Fourier 
phase (j>2i places this Cepheid away from the fundamental 
mode sequenc e and possibly onto th e first overtone sequence. 
On this basis, iKienzle et al.l l|l999l) proposed that V440 Per 
is an overtone pulsator. This hypothesis was further sup- 
ported by determination of the phase-lag between the light 
curve and the radial velocity curve, A"l>i = cj>Y^ — 0™"®, 
which pl aced V440 Per away from the fundamental sequence 
as well l|Ogloza et al.ll2000l ). However, V440 Per is a smaU 
amplitude, nearly sinusoidal variable. Even the b est then 
available radial velocity data (|Burki fc Benj[l982l ) yielded 
large errors of the 1st harmonics Fourier parameters. The 
measurement error of the phase-lag A$i was large, too. 
This, and i nconsistency with their hydrodynamic pulsation 
models led'Sz abo et al.l (|2007l ) to dispute the mode identifi- 
cation of V440 Per. They noted that the membership of this 
Cepheid in the fundamental mode sequence could not be re- 
jected at the 3(T confidence level. They argued that V440 Per 
is not an overtone pulsator, but rather a fundamental mode 
Cepheid of a very low amplitude. 

In the present paper we report results of an extensive 
campaign of the observations of V440 Per with the Poznan 
Spectroscopic Telescope, lasting almost a year and aiming 
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Table 1. PST radial velocities of V440 Per. 
Full table in the electronic form only. 

MJD A?;"' ^v'^.corr 



54327.0490 


-0.61 


-0.52 


54332.0517 


3.04 


3.14 


54332.0683 


2.96 


3.06 


54649.0009 


2.28 


2.12 


54650.9302 


2.75 


2.57 


54650.9419 


2.76 


2.58 



a) zero shifted by about —26.3 ± 0.2 km/s 

b) corrected for instrumental effects (see text) 



at obtaining high quaUty data suitable for the detailed diag- 
nostic and comparison with the nonlinear models. In Sect .[2] 
we describe our instrument and observations. In Sect.Owe 
discuss quality of our radial velocities and derive Fourier pa- 
rameters of V440 Per. Comparison with Fourier parameters 
of other Galactic Cepheids and identification of the pulsa- 
tion mode of V440 Per is discussed in Sect.|4] In Sect.[5]we 
perform detailed comparison of V440 Per with the nonlinear 
overtone Cepheid models. Our conclusions are summarized 
in Sect.ini 



2 OBSERVATIONS AND DATA REDUCTION 

Our observations were obtained with the new Poznari Spec- 
troscopic Telescope (PST) of Adam Mickiewicz University. 
Its full description is not published yet, hence we devote 
some space here to instrument description and data quality 
evaluation. 



2.1 Poznaii Spectroscopic Telescope 

PST is located in Poland at Borowiec station, 20 km south 
from Poznari city, at a meagre elevation of 123 m above see 
level. The PST consists of parallel twin 0.4 m Newton tele- 
scopes of the f-ratio 4.5, fixed on a single parallactic fork 
mount. An acquisition box at Newton focus of each tele- 
scope holds the tip of the fiber feeding our spectrograph, 
the Thor/ Argon calibration lamp and the autoguider cam- 
era SBIG ST-7. 

The telescope feeds via a fiber an Echelle spectrog raph, 
a clone of MUSICOS design l|Baudrand fc Bohm|[l993 ). red 
arm only. Our spectra are recorded with the low noise An- 
dor DZ436 camera fitted with 2fcx2fc E2V 42-40 back illumi- 
nated CCD chip, cooled with Peltier cells. About 60 orders 
are recorded, covering a spectral range of 4480 -9250 A at 
the inverse resolution of A/(5A = 35000. The spectrograph is 
located in a thermo-isolated enclosure in the telescope dome. 
The sliding-roof dome, the telescope and its spectrograph all 
operate under full computer control. A full description of our 
instrumentation will be published elsewhere. Our system is 
operated interactively from the terminal. Such an operation 
mode does not compromise our data quality but requires 
excessive work-load. A substantial software effort is needed 
to achieve the fully robotic operation. 
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Figure 1. Radial velocities of V440 Per vs. phase. First and sec- 
ond order Fourier fits are plotted with dashed and solid lines, 
respectively. The error bar indicates ± 1 standard deviation of 
residuals. The velocities are corrected for small instrumental ef- 
fects (see text). Zero point is arbitrary. 
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Figure 2. Top : residuals of the first order Fourier fit of V440 Per 
radial velocities (dashed line in Fig.[TJ vs. phase. 1st harmonic 
is very small {A2 = 140 ± 15m/s), yet clearly visible. Bottom: 
residuals of the second order Fourier fit vs. time. 



2.2 Data pipeline 

Routine CCD reductions up to spectra extraction, wave- 
length cafibration and velocity measurement are performed 
using IRAF tasks combined into our reduction pipeline. Ve- 
locities are measured by cross-correlation, using IRAF FX- 
COR task. The internal error estimates from FXCOR, of 
order of few km/s, relate to the line width and not to the 
actual measurement precision. They may serve for weight- 
ing purposes, though. So far, we employed no standard star 
calibrations and our velocities are measured solely with re- 
spect to the Th/Ar lamp. For short span observations such 
a primitive procedure still yields root-mean-square (RMS) 
residuals in the 100-200m/s range, for stars brighter than 
11 mag. As demonstrated in Sect.|3l for intensively observed. 
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Table 2. Fourier parameters of V440 Per radial velocity curve 



Name 


Value 


Error 


Unit 


To 


54498.783 


0.008 


MJD 


P 


7.5721 


0.0006 


d 


Ai 


2.480 


0.015 


km/s 


A2 


0.140 


0.015 


km/s 


R21 


0.056 


0.006 


1 


4>21 


2.759 


0.117 


rad 



strictly periodic stars, any long-term effects of a floating in- 
strumental zero-point may be removed with some assurance, 
yielding again RMS residuals less than 150 m/s, over one 
year span of data. 



2.3 Observations 

In total 158 radial velocity measurements of V440 Per were 
obtained with PST from 2007 August 15 till 2008 July 03. To 
reach signal-to-noise ratio of 70 we exposed spectra for 10 - 
15min and up to several spectra were obtained per night. 
The observed velocities are listed in Table [T] published in the 
electronic form. Phase-folded data with fitted Fourier series 
(Sect. [3]) are displayed in Fig.[T] In Fig.[2]we plot residuals of 
the first order and the second order fits vs. pulsation phase 
and vs. time. Inspection of the plots demonstrates that our 
coverage was reasonably uniform, both in time and in fre- 
quency. From the residuals we estimate standard error of 
our individual measurements as 130 m/s. 



3 FOURIER PARAMETERS 

The plot of the phase-folded radial velocities reveals near si- 
nusoidal variations with peak-to-peak amplitude of 5 km/s. 
Some data consistency checks are due prior to drawing any 
final conclusions. Cepheid phases and frequencies are known 
to vary slightly. Additionally, the instrument stability over 
9 months needs checking, too. A preliminary nonlinear least 
squares fit of our data with the Fourier series of three har- 
monics terms (third order fit) yielded no significant 2nd 
harmonics of the main frequency nor the period derivative 
term. Our fitted frequency of 0.13206 ± 0.00001 c/d is less 
accurate, yet consistent within errors with the frequency 
derived by combining our dat a with the earli e r mea sure- 
ments of iBurki fc Benz (l98^, lArellano Feirol (|l984l ') and 
iGorvnva et al.l l|l992. 1996, 199^). Comparison of the data 
shows, that the weighted zero point shift of our velocities 
with respect to the previous authors is —26.3 ± 0.2 km/s. 

The RMS deviation from the Fourier fit of all our data 
was 164 m/s, well in excess of 135 m/s obtained for the first 
part of the dataset. Worse, inspection of the residuals plot- 
ted against time sometimes revealed non-gaussian, bi-modal 
distribution. The origin of both effects seems to be instru- 
mental. To confirm that, we expanded our Fourier model by 
including two linear terms proportional to the time inter- 
val from the mid-epoch and to the hour angle of the star 
at the moment of observation. The fitted values of these in- 
strumental correction coefficients were significant at 6a and 
5a level, respectively. From the overall covariance matrix we 
find maximum absolute value of the correlation coefficients 
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Figure 3. Fourier phase (/)2i vs. pulsation period for Cepheid ra- 
dial velocity curves. Fundamental mode Cepheids displayed with 
asterisks, except for low amplitude ones (Ai < 10 km/s) marked 
with open triangles. Overtone Cepheids dis played with fil l ed cir - 
cles. Observational data p oints taken from iKienzle et al.l (Il999l ) 
and iMoskalik et al.l l|2009h . Straight line indicates the best fit to 
the segment of the fundamental mode progression. V440 Per dis- 
played with open circle. 



of 0.30, consistent with little interference between different 
fitted terms. These instrumental corrections reached up to 
±100 m/s. At this stage, we have no explanation for these 
corrections. 

Our final second order Fourier fit, supplemented with 
the instrumental correction terms, yielded RMS deviation 
of 130 m/s, consistent with that obtained from the short- 
span observations. The values of the Fourier parameters of 
V440 Per radial velocity curve are listed in Table [5] (for 
exact formulae defining Fourier parameters and their errors 
see Appendix [A}. 



4 PULSATION MODE OF V440 PER 

Pulsation mode of a Cepheid can be esta blished by mea- 
surin g Fourier phase </)2i of its li ght curve llAntonello et al.l 
Il99(l ) or radial velocity curve (jKienzle et all ligggl ). This 
Fourier parameter does not depend on the pulsation am- 
plitude of the star and for each mode it follows a different, 
tightly defined progression with the pulsation period. For pe- 
riods at which the two 021 progressions are well separated, 
secure mode identification can be achieved. 

In Fig.|3] we plot velocity <j>2i of short period Galac- 
tic Cepheids against their pulsation period P. Fundamen- 
tal mode pulsators and overtone pulsators are displayed 
with different symbols. V440 Per is plotted with an open 
circle. It is immediately obvious, that it is located far 
apart from the fundamental mode progression. This notion 
can be put on a quantitative basis. We selected a sample 
of 23 fundamental mode Cepheids with periods P in the 
range of Po ± 1.1 day, where Po = 7.5721 day is the pe- 
riod of V440 Per. To this sample we fitted a straight line 
021 (f) = a,{P — Po) + b. With this procedure, we find that 
at the period of V440 Per the expected 02i of the funda- 
mental mode Cepheid is 3.628 rad, with the average scatter 
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Table 3. Convective parameters of pulsation models discussed in 
this paper, as, «c, ctp a^nd 7r are given in the units of standard 
values (see Smolec & Moskalik 2008). In the last column we give 
a linear upper limit for first overtone period (Pmax). 



12 16 
Aj [km/s] 

Figure 4. Residuals from the mean fundamental mode P — 4>2i 
relation vs. pulsation amplitud e A i. Observationa l data points 
taken from lKienzle et al.l Il999l) and lMoskalik et al.lll2009^ . Same 
symbols as in Fig.|3] 



of individual values of ao = 0.026 rad. This estimate of the 
intrinsic scatter is conservative, as the nominal 021 mea- 
surement errors would account for at least half of it. The 
021 value measured for V440 Per is 2.759 rad, with an error 
of (71 = 0.117 rad. The probability distribution of the 02i 
offset. A, is obtained by convolution of two normal distribu- 
tions, N{0,ao) and N{A,ai). It may be demonstrated, that 
the result is another normal distribution A'^(A, \J (y^-\- af ), 
essentially by virtue of the law of error combination. By sub- 
stituting A = 3.628 - 2.759 = 0.869 rad, we find that the 
observed velocity 021 of V440 Per deviates from the fun- 
damental mode sequence by 7.25a. Thus, from the purely 
observational evidence we conclude that V440 Per does not 
pulsate in the fundamental mode. Consequently, it must be 
an overtone Cepheid. 

The argument presented above critically depends on the 
assumption that velocity 02i (at a given period) does not 
depend on the pulsation amplitude, which for V440 Per is 
lower than for the fundamental mode Cepheids used as com- 
parison. Such an assumption is well justified by numerical 
computations (i.e. Buchler et al. 1990; Smolec & Moska- 
lik 2008), nevertheless it should be verified with available 
data. In Fig.|3]we plot residuals from the mean fundamental 
mode P — 021 relation vs. amplitude Ai. In this case the 
mean relation is defined as the parabola fitted to all fun- 
damental mode Cepheids displayed in Fig. [31 For Ai in the 
range of 10 — 17km/s there is no correlation of 02i residu- 
als with pulsation amplitude. For Ai < 10 km/s there is a 
weak evidence of 02i increase. As the data are very scarce, 
we are not convinced that this increase is significant. If it 
were real, however, it would only strengthen our conclusion 
that V440 Per strongly deviates from the fundamental mode 
sequence. 



5 MODELING OF V440 PER 

Classification of V440 Per as a first overtone pulsator, dis- 
cussed in the previous Section, is based solely on the mor- 
phological properties of Cepheid velocity curves. This mode 
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identification can be strengthened by comparing the ve- 
locity curve of V440 Per with those of hydrodynamical 
Cepheid mod els. Such a comparison was first performed by 
iKienzle et al.l (11999), who used the unpublished radiative 
models of Schaller & Buchler (1994). They showed that the 
theoretical progression of velocity 02i supported the over- 
tone classification of V440 Per. However, their conclusion 
was somewhat weakened by a large error of then available 
021 of this star. 

In the present section, we confirm and extend the re- 
sults of iKienzle et al.l (|l999l ). It is easy to show, that the 
velocity curve of V440 Per is incompatible with the funda- 
mental mode Cepheid models. Indeed, all published models 
display velocity 02i higher that 3.0 rad at all periods and 
higher than 3.5 rad at the period of V440 Per {i.e. Buchler 
et al. 1990; Moskalik et al. 1992; Smolec & Moskalik 2008). 
This holds true both for the convective models and for the 
older radiative models. The computed velocity 02i is very 
robust and shows no sensitivity to the treatment of convec- 
tion, choice of opacities or to details of the numerical code. 
Most importantly, it is insensitive to the pulsation amplitude 
(c/. Figs. 8, 11 of Smolec & Moskalik 2008). Clearly, the only 
chance to match the observed velocity curve of V440 Per is 
to search for an appropriate overtone model. 

With this goal in mind, we computed several sequences 
of the convective overtone Cepheid models. We show that 
V440 Per fits theoretical 02i progression of the first overtone 
Cepheids and that its velocity Fourier parameters can be 
accurately reproduced with hydrodynamical computations. 

Modeling of such a long period overtone pulsator is 
not an easy task. The satisfactory models have to repro- 
duce Fourier parameters and the long period of this vari- 
able. The current hydrocodes used for modeling of the radial 
pulsations adopt a time-dependent convection models (e.g. 
Stellingwerf 1982, Kuhfufi 1986). These models introduce 
several dimensionless parameters, that should be adjusted 
to match the observational constraints. V440 Per with its 
exceptionally long period, offers an opportunity to impose 
interesting limits on the free parameters of the convection 
treatment. 

In our computa t ions w e use the convecti ve hydro c ode o f 
ISmolec fc MoskalikI (|2008D . The code adopts iKuhfuS (|l986D 
dynamical convection model reformulated for use in stellar 
pulsation calculations. Kuhfufi model is physically well mo- 
tivated and self-consistent. It contains 8 scaling parameters, 
summary of which, together with the detailed description 
of the model equations is provided by ISmolec fc MoskalikI 
(|2008l ). The values of parameters used in the present paper 
are given in Table [3] All our Cepheid models are constructed 
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Figure 5. Theoretical instability strips for convective parameters 
of set A. 



in the way described in ISmolec fc MoskalikI (|2008l ). In all 
computations we use Galactic chemical c omposition (X = 
0.70, Z = 0.02) and O PAL opacities Jldesias fc Rogerl 
1 19961 ) computed with the lGrevesse fc Noelsl | T993[) mixture. 
We use the m a ss-lum inosity (M — L) relation derived from 
ISchaller et al.l l| 19921 ') evolutionary tracks (log(I//Z/0) = 
3.56 log(M/M0) +0.79). 



5.1 Pulsation period of V440 Per 

Theoretical instability strips (IS) computed with the linear 
pulsation code for set A of the convective parameters are 
presented in Fig.[S] Instability strips for the first overtone 
and for the fundamental mode are enclosed with the thick 
and thin fines, respectively. The dotted lines correspond to 
constant values of the first overtone period, as indicated in 
the Figure. 

First overtone IS does not extend toward an arbitrar- 
ily high luminosity. The linear computations yield an up- 
per limit for the first overtone pulsation period, reach- 
ing ~ 8.1 day in Fig. [5] (see also Table [Sjl. This limit de- 
pends mostly on the adopted convective parameters, but 
also on the M — L relation and on the metalicity. However, 
to find a model, that satisfies V440 Per period constraint, 
it is not enough to assure appropriate linear period limit 
(P > 7.5 day). This is because the maximum overtone pe- 
riod at full amplitude pulsation is determined by the nonlin- 
ear effects. 

Between the linear blue edges of the first overtone and 
of the fundamental mode and below their crossing point, the 
first overtone pulsation is the only possibility. As one can see 
in Fig. [5] such a region is very narrow. A significant part of 
the overtone IS lays inside the fundamental mode IS. For 
higher luminosity and overtone periods longer than 6 days, 
the instability strips of both modes entirely overlap. 

In the region, where both modes are linearly unstable. 



the final pulsation state is determined by the nonlinear ef- 
fects (modal selection). Full amplitude pulsation in one of 
the modes usually suffices to saturate the instability. In prin- 
ciple, a double-mode pulsation is also possible, but it is very 
unlikely at long overtone periods that are of interest in this 
paper (no double-mode pulsator with such a long overtone 
period is known either in the Galaxy, or in the Magellanic 
Clouds; see Soszyriski et al. 2008, Moskalik fc Kolaczkowski, 
2008). Therefore, construction of models matching period of 
V440 Per requires both linear and nonlinear computations 
and can be used to constrain parameters of convection model 
we use. 

It is hard to const rain all 8 parame ters entering the 
model. As was shown bv lYecko et al.l l| 19981 ). different combi- 
nations of the convective a parameters may yield essentially 
the same results. Therefore, we decided to freeze four param- 
eters, for which standard values are in use. In all studied sets 
of the convective parameters (Table |3]), mixing- length is set 
to a = 1.5 and for Oa, etc and ad we use the standard values 
(see Smolec & Moskalik 2008 for details). Set A represents 
the simplest convection model, without turbulent pressure, 
turbulent fiux or radiative losses. These effects are turned 
on in set B (radiative losses) and in set C (turbulent pres- 
sure and turbulent fiux). One of the crucial factors of the 
convection model is eddy-viscous damping, the strength of 
which is determined by the Ofm parameter. The lower the 
eddy- viscous damping is (lower am), the more linearly un- 
stable become the models and, consequently, the higher are 
their pulsation amplitudes. 

The most interesting outcome of the linear computa- 
tions are the period limits, Pmax, given in the last column of 
Table |3] At the linear theory level, all sets of convective pa- 
rameters listed in the Table can satisfy V440 Per period con- 
straint. This was assured by adjusting eddy- viscous damping 
parameter, am- In case of sets A and C (no radiative losses) 
low values of am are required. 

For each parameter set we computed a sequence of non- 
linear models, running parallel to the first overtone blue 
edge, at a constant distance AT. In this Section we use 
AT = 75 K. This choice is arbitrary, but we believe, that such 
models should reproduce most of the observed first over- 
tone variables. The static models were initialized in the first 
overtone mode (as described in Smolec fc Moskalik 2008) 
and their time evolution was followed, until the final pulsa- 
tion state was reached. The computed radial velocity curves 
were then decomposed into Fourier series. In Fig. [S] we plot 
amplitude A\ vs. pulsation period for three sequences of 
models. The most satisfactory results are obtained for set B 
of the convective parameters. The pulsation amplitudes for 
this sequence are moderate and they match well the obser- 
vations (this conclusion will be strengthened in Sect. 15. 2| ). 
Also, for this sequence we obtain the overtone models with 
the longest periods. Sequences computed with the param- 
eter sets A and C extend to periods not exceeding 6 days. 
All more massive models {i.e. with longer overtone periods) 
switched into the fundamental mode pulsation and therefore 
are not plotted in the Figure. The amplitudes of the mod- 
els computed with the parameter sets A and C are larger. 
This is due to the adopted low value of the eddy-viscous 
damping, needed to obtain high enough limit for the linear 
overtone period. Any further decrease of am does not help. 
Although Pmax increases, yet models switch into the funda- 
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Figure 6. Fourier amplitude A\ for three sequences of overtone 
Cepheid models. Model amplitudes scaled by constant projection 
factor of 1. 4. Except of V440 P er, o bservational data poin ts are 
taken from lxienzle et al.l lll999h and lMoskalik et al.l ||2009|) . 



mental mode at periods well below of 7.57 days required to 
match V440 Per. Additionally, at the short periods the am- 
plitudes of the first overtone models become unacceptably 
high. We conclude, that inclusion of the radiative losses is 
necessary to reproduce the long overtone period of V440 Per. 
As we will demonstrate in Sect . 15.21 with the set B of the 
convective parameters we can also reproduce the Fourier pa- 
rameters of V440 Per, as well as the overall progression of 
Fourier parameters for all observed overtone Cepheids. 

Our co n clusio n is consistent with the results of 
ISzab6 et all ()2007l ). For convection with the radiative 
losse^, their models extend toward high overtone periods 
(parameter set A of Szabo et al. 2007). Without radiative 
losses (set B of Szabo et al. 2007) and with much higher 
eddy-viscous dissipation than in our set A, their longest 
overtone periods fall below 3 days. Co mparing our mod- 
els wi th the overtone Cepheid mo dels of iFeuchtinger et all 
(|2000l ) we note some inconsistency. IFeuchtinger et al.l (|200(]f ) 
obtain long overtone periods for all sets of convective param- 
eters considered in t heir paper. Their results are inconsistent 
with the results of lSzabo et al.l |2007j) (set A of Feuchtinger 
et al. 2000 and set B of Szabo et al. 2007 are identical, but 
their linear results are significantly different) as well as with 
our results. We trace this discrepancy to different evalua- 
tion of the superadiabatic gradi ent, y = V — Vg , in ou r 
and in the Vienna code used by IFeuchtinger et al.l (|2000h . 
In the Vienna code, the radiative pressure contribu tion is 
neglected in computation of Va (|Feuchtingetlll999l l. This 
leads to higher values of Va {~ 0.4 far from partial ioniza- 
tion regions, see Figs. 3 and 4 of Wuchterl & Feuchtinger 
1998) and lower values of the superadiabatic gradient. Con- 
sequently, in their models convection is less pronounced and 
the models are more unstable. In fact, we were able to repro- 
duce their linear results when we neglected the contribution 
of radiation in the evaluation of Va. 



5.2 Fourier parameters of V440 Per 

In this Section we focus our attention on set B of the convec- 
tive parameters. As demonstrated in the previous Section, 



^ Note how e ver, t hat in Florida-Budapest hydrocode used by 
ISzabo et al radiative losses are modeled in a different 

way (see also Smolec & Moskalik 2008). 
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Figure 7. Fourier parameters of Cepheid radial velocity curves. 
Model sequences are computed at constant distances from the 
first overtone blue edge, as indicated in the Figu re. Except of 
V440 Per, observational data poi nts are taken from lKienzle et al.l 
(1999!) and lMoskalik et al.l ('2009). Error bars for most of the am- 
plitudes are smaller than the symbol size and are not plotted. 



this parameter set allows to reproduce the long overtone pe- 
riods and the observed pulsation amplitudes. In Fig. [7] we 
compare the velocity Fourier parameters of Galactic over- 
tone Cepheids with those of the hydrodynamical models 
computed with the parameters of set B. Four sequences of 
models displayed in the plot run at different distances from 
the first overtone blue edge (AT = 25 K, 75 K, 125 K, 175 K). 
The overall observed progression of Fourier parameters is 
well reproduced with our models. For 02i, we notice that 
the model sequences are shifted toward somewhat shorter 
periods. This can be easily explained. The characteristic pro- 
gression of velocity 021 with period is caused by the 2:1 res- 
onance between the first overtone and the fourth overtone 
llKienzle et al. I I1999I : IFeuchtinger et al.ll2000D . The center of 
this resonance is located at P = 4.2 — 4.6 day. The exact 
location of the resonance in hydrodynamic models depends 
mostly on the chosen mass-luminosity relation. This relation 
was not adjusted in our calculations to match the observed 
resonance progressior(f|. As shown in Fig.O the resonance 
center (displayed with the dashed line) crosses the center of 
the first overtone IS for periods shorter than 4 days, which 
explains the horizontal shift in Fig. [7] Taking into account 
this shift, amplitudes, amplitude ratios and Fourier phases 
agree satisfactorily with observations, although at the short 



^ Adjustment of M — L relation to satisfy the resonance con- 
straint is not an easy task. Both the slope and the zero point 
of the adopted M—L relation affect location of the resonance and 
the shape of the 4'2l progression. In a forthcoming publication 
(Smolec & Moskalik, in preparation) we address this problem in 
detail. 
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periods theoretical (f)2i values are slightly too high. Consid- 
ering V440 Per, we see the very good agreement with the 
model sequence closest to the blue edge of the overtone in- 
stability strip {AT — 25 K). An exact match can be easily 
obtained for the model sequence located slightly closer to 
the blue edge. We note, that the progression of the overtone 
(j)2i at long periods is insensitive to the choice of AT. The 
models predict slightly higher 021 than actually observed in 
V440 Per, but the two values are consistent within the error 
bax. 
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6 CONCLUSIONS 

Using our new Poznah Spectroscopic Telescope we obtained 
158 high precision radial velocity observations of a low am- 
plitude Cepheid V440 Per. We constructed the pulsation 
velocity curve of V440 Per and we were able to reliably de- 
tect its 1st harmonic, with amplitude of only 140 ± 15m/s. 
The measured Fourier phase (f)2i = 2.76 ± 0.12rad differs 
from the values observed in fundamental mode Cepheids 
of a similar period by 7.25a". Thus, we demonstrated on 
purely morphological ground that V440 Per does not pulsate 
in the funda mental mo de. This settles the d ispute between 
ISzabo et all ([2007) and iKienzle et aD ljl999h and allows to 
classify V440 Per as an overtone pulsator, the one with the 
longest period identified so far {P = 7.57 day). 

Our results demonstrate that with suitable care, our in- 
expensive instrument featuring MUSICOS Echelle spectro- 
graph and a small robotic telescope can achieve stability and 
precision surpassed only in the extrasolar planet searches. 
Note, that we employed neither the iodine cell nor environ- 
ment control. Yet, our observations prove that the secure 
mode identification is feasible even for the very low ampli- 
tude Galactic Cepheids. 

The overtone pulsation of V440 Per has interest- 
ing theoretical consequences. To investigate them we em- 
ployed our convectiye line ar and nonlinear pulsation codes 
(|Smolec fc Mo skalik 2008). The first overtone linear models 
are already constrained by the value of the pulsation period 
at the center of the 2uji = toi resonance (jKienzle et al.ll 19991 : 
iFeuchtinger et al1l2000l ') . The very existence of V440 Per im- 
poses additional constraints. Namely, to linearly excite the 
first overtone and then to obtain a nonlinear full-amplitude 
overtone pulsation of such a long period, one has to fine 
tune the dynamical convection model used in the pulsation 
calculations. Our numerical experiments demonstrate that 
in the convective energy transport radiative losses must be 
properly accounted for to maintain consistency with the ob- 
servations. With this effect taken into account, the nonlinear 
overtone Cepheid models not only reproduce the exception- 
ally long period of V440 Per, but they also reproduce neatly 
all Fourier parameters of its pulsation velocity curve. No 
such agreement can be achieved with models pulsating in 
the fundamental mode. These results of the hydrodynam- 
ical modeling provide additional support for our empirical 
classification of V440 Per as a first overtone pulsator. 
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APPENDIX A: FOURIER COEFFICIENTS AND 
THEIR ERRORS 

It seems useful to collect in one place all the relevant formu- 
lae for Fourier coefficients of velocity v{t) and their errors 
a. Defining Fourier series in the following way 

iV 

v{t) - ao = An sm{nujt + 

n = l 
JV 

= ^ ^ (on cos nut + bn sin nmi) (Al) 

n = l 

where A'^ is the order of the fit, we obtain 

/ A 

An = Vffln + bl and Rni = -r^ (A2) 

A\ 

tan<^„ = ^ and cj>ni = (pn ~ n(j>i (A3) 

On 

From variations of these equations we find 



SAn = J (A4) 



a„5a„ + bnSbn aiSai+biSbi 



SRnl = -Rnl ( T2 J2 ) (A5) 



6,1 (5a,, — a„56„ 



(A6) 



fe„5a„ - a„56„ 6i5ai— ai56i 

d^nl = 72 42 (A7) 

where RHS are scalar products of the vector 
((5a„, 5ai, (5bi) = and the corresponding deriva- 
tives {d / dan,d / dbn,d / da\,d / dbi) = L[-]. For example 
S(j>ni ~ L[(j)ni] ■ SF and third component of vector I/[(;/)„i] 
is d(j)ni/dai — —nbi/Al. By the law of error propagation, 
variance of (jini is 

(T^^j = Var[4>„i] = L[(l)„i]xCov[an,b„,ai,bi]xL[4>„i]'^{A8) 

where Cov[-] denotes the covariance matrix of raw Fourier 
coefficients (Ea JAlf) . x denotes matrix multiplication and 
superscript T indicates matrix transpose. Analogous expres- 
sions hold for other coefficients. The complete linearized 
transformation matrix A4 is obtained by substitution of the 
L vectors corresponding to Eqs. HA4HA7|) as its rows. The 
full covariance matrix of all Fourier coefficients is seldom 
needed in practice. It may be obtained from Eq. (|A8|I by 
substitution of M for L. 



